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A  strongly  disturbed  layer  of  ionization  irregularities  that  is 
used  as  a  propagation  channel  for  radio  waves  can  degrade  the  propagating 
wave  and  thereby  affect  the  resulting  measurements  at  the  receiving 
antenna.  The  antenna  aperture  itself  also  affects  measurements  of  the 
received  signal  Sy  its  inherent  averaging  process.  Here  an  analytic  solu¬ 
tion  for  the  two-position,  two-frequency  mutual  coherence  function,  valid 
in  the  strong-scatter  limit,  is  used  to  characterize  the  propagation 
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20.  Abstract  (Continued) 

channel,  lhe  channel  itself  consists  of  a  thick  slab  of  anisotropic  elec¬ 
tron  density  irregularities  that  are  elongated  in  the  direction  parallel 
to  the  earth's  magnetic  field. 

Analytic  expressions  are  obtained  that  give  the  effect  of  the 
aperture  antenna  on  measurements  of  received  power,  decorrelation  time  (or 
distance),  mean  time  delay,  time  delay  jitter  and  coherence  bandwidth. 
These  quantities  are  determined  as  functions  of  the  aperture  diameter  and 
of  the  angle  between,  the  magnetic  field  and  the  direction  of  propagation. 
It  is  shown  that  in  strong  turbulence  aperture  averaging  can  be  a  signi¬ 
ficant  factor  in  reducing  the  received  power  by  angular  scattering  loss, 
increasing  the  observed  signal  decorrelation  time  via  aperture  averaging, 
and  reducing  the  time  delay  jitter  by  suppression  of  signals  received  at 
off-boresight  angles. 

Results  are  presented  for  two  cases.  One-way  propagation 
through  an  ionospheric  communication  channel  is  considered  where  both 
transmitter  and  receiver  utilize  aperture  antennas.  This  result  is  easily 
extended  to  the  case  that  one  of  the  antennas  is  an  omnidirectional  point 
source,  corresponding  to  the  usual  case  of  trans ionospheric  satellite 
communication  from  a  small  satellite  antenna  to  a  large  ground  based 
receiver.  The  second  case  involves  transmission  and  reception  of  a  radar 
signal  that  travels  through  a  disturbed  ionospheric  channel  to  a  target 
located  in  free  space.  This  case  is  applicable' to  the  situation  of  a 
large  antenna  aboard  a  space  based  radar  or  to  the  case  of  a  ground  based 
defense  radar. 
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SUMMARY 


In  this  report  analytic  expressions  are  derived  for  the  effect 
of  aperture  antennas  on  measurements  of  signals  that  have  propagated 
through  strong  anisotropic  turbulence.  Results  are  given  for  the  effect 
of  Gaussian  apertures  on  measurements  of  received  signal  power,  decorrela¬ 
tion  distance  (or  time),  mean  time  delay  and  time  delay  jitter.  The 
geometries  considered  correspond  to  the  case  of  a  single,  one-way  propaga¬ 
tion  path  between  two  aperture  antennas,  located  in  free  space  and 
separated  by  a  layer  of  turbulence  artd  to  the  case  of  a  monostatic  radar 
where  radar  and  target  are  on  opposite  sides  of  a  strong  scattering  layer. 

All  calculations  depend  on  the  analytic  solution  for  the  two- 
position,  two-frequency  mutual  coherence  function  for  spherical  wave 
propagation  in  an  anisotropic  layer  of  electron  density  irregularities. 

The  analytic  solution  for  propagation  in  a  thick  layer  is  derived  here 
using  the  quadratic  phase  structure-function  approximation  valid  for 
strong  turbulence.  This  result  is  then  specialized  to  a  thin  phase-screen 
approximation  to  facilitate  analytic  calculation  of  the  effects  of 
apertures. 

It  is  shown  that  antennas  that  are  larger  in  diameter  than  the 
decorrelation  distance  that  would  be  measured  by  an  omnidirectional 
antenna  can  experience  significant  angular  scattering  loss  and  exhibit 
increased  measurements  of  decorrelation  distance,  and  decreased  measure¬ 
ments  of  mean  time  delay  and  time  delay  jitter.  Increased  measurements  of 
signal  decorrelation  distance  are  caused  by  the  averaging  effect  of  the 
antenna  aperture  that  smooths  (eliminates)  the  small  scale  signal  fluctua- 


tions.  Decreased  measurements  of  mean  time  delay  and  time  delay  jitter 
are  caused  by  the  action  of  the  aperture  to  cut  off  signal  contributions 
that  originate  from  off-boresight  angles  and  have  therefore  experienced 
more  time  delay  than  experienced  over  more  direct  signal  paths. 

Generally,  the  effects  of  aperture  averaging  first  begin  to  be 
significant  for  aperture  sizes  .that  approach  ten  times  the  decorrelation 
distance  measured  with  a  point  antenna.  For  even  larger  apertures,  the 
aperture  has  a  significant  effect  on  the  received  signal  at  the  antenna 
output  that  can  be  computed  using  the  formulas  and  curves  provided  herein. 
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SECTION  1 
INTRODUCTION 


Large  .'nigh-gain  antennas  are  used  in  many  radar  and  communica¬ 
tions  applications  to  increase  the  energy  collected,  to  increase  angular 
accuracy,  and  to  provide  protection  against  jamming.  If  the  wavefront  at 
the  antenna  aperture  experiences  scintillation,  large  apertures  can  act  to 
average  the  signal  and  thereby  modify  the  observed  signal  properties.  In 
the  late  fifties  Wheelon  (1957,  1959)  obtained  the  well-known  aperture 
smoothing  effect  for  the  measurement  of  phase  fluctuations  by  a  finite 
circular  aperture.  Since  then,  various  investigators  have  studied  the 
aperture  smoothing  effect  on  intensity  scinti 1  lotions  in  the  optics  regime 
(Fried,  1967;  Tatarskii,  1971;  Homstad  et  al.,  1974)  and  also  the  related 
effect  of  optical  beam  wave  oropagation  to  an  infinitesimal  receiver  (Lee 
and  Harp,  1969;  Ishimaru,  1969).  Knepp  (1975)  obtained  results  for  the 
effects  of  an  antenna  aperture  on  measurements  of  the  in-phase  and  quad¬ 
rature  components  and  their  spatial  derivatives,  for  the  case  of  weak 
scattering. 

Scintillation  or  rapid  variations  in  the  amplitude,  phase,  and 
anglc-of-arri val  of  a  propagating  wave  is  often. observed  over  satellite 
links  at  VHF  and  UHF .  Scintillation  can  also  be  observed  at  frequencies 
as  high  as  the  GH^  range  (Pope  and  Frit?,  1971;  Skinner  et  al.,  1971; 

Taur,  1976)  and  is  occasionally  severe,  even  at  L-band  (Fremouw  et  al., 
1978).  However,  aperture  averaging  is  not  generally  important  at 
satellite  frepiencies  (UHF -GHz )  under  ambient  or  naturally  perturbed 
ionospheric  conditions.  - 


Worst  case  or  Rayleigh  amplitude  scintillation  is  likely  to 
occur  if  the  ionosphere  is  highly  disturbed,  as  for  example  by  high  alti¬ 
tude  nuclear  explosions  (Arendt  and  Soicher,  1964;  King  and  Fleming,  1980) 
or  by  chemical  releases  (Davis  et  al.,  1974;  Wolcott  et  al.,  1978). 
Increased  electron  concentrations  and  the  irregular  structure  of  the  ioni¬ 
zation  can  lead  to  intense  Rayleigh  signal  scintillation  at  frequencies  as 
high  as  the  7-8  GHz  SHF  band  (Knepp,  1977).  Consequently,  the  effects  of 
scintillation  are  important  to  any  UHF  through  SHF  communications  or  radar 
system  that  must  operate  through  an  ionospheric  channel  and  that  may  have 
to  operate  in  highly  disturbed  propagation  environments. 

For  cases  of  severe  scintillation  where  the  signal  varies  across 
the  area  of  the  receiving  aperture,  the  effect  of  the  aperture  can  be 
significant.  It  is  well  known  that  an  aperture  antenna  acts  to  cut  off 
energy  that  is  incident  at  off-boresi ght  angles  where  the  antenna  gain  is 
reduced.  This  effect  may  also  be  viewed  as  the  result  of  averaging  or 
coherent  processing  of  the  electromagnetic  field  intercepted  at  the 
aperture  location.  In  this  report  th,e  effects  of  aperture  averaging  are 
analytically  calculated  for  the  case  where  Gaussian  antenna  beams  are  used 
by  transmitter  and  receiver.  Results  are  presented  for  two  different 
physical  situations.  One-way  propagation  through  an  ionospheric  communi¬ 
cation  channel  is  considered. where  both  transmitter  and  receiver  utilize 
aperture  antennas.  This  result  is  easily  extended  to  the  case  that  one  of 

- - the-antennas  is  a  point  source,  corresponding  to  the  usual  case  of  trans- 

ionospheric  satellite  communications.  The  second  antenna  geometry  corre¬ 
sponds  to  transmission  and  reception  of  a  radar  signal  that  travels 
through  an  ionospheric  channel  to  a  target  located  in  free  space.  This 
corresponds  to  the  situation  of  a  long-range  ground  or  space  based  radar. 
It  is  shown  that  aperture  averaging  can  greatly  affect  measurements  of 
scattering  loss,  signal  decorrelation  time,  mean  time  delay  and  time  delay 
jitter.  Simple  analytic  expressions  are  given  for  all  these  quantities  in 
terms  of  the  geometry  of  the  propagation  path  and  the  severity  and  struc¬ 
ture  of  the  ionization  irregularities. 
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To  obtain  results  applicable  to  a  general  geometry  of  the  line- 
of-sight  relative  to  the  field  aligned  ionization  structure,  it  is  neces¬ 
sary  to  obtain  an  analytic  solution  for  the  two-position,  two-frequency, 
mutual  coherence  function  (MCF)  for  spherical  wave  propagation  through  a 
thick  layer  of  anisotropic  electron  density  irregularities.  It  is  assumed 
here  that  strong  scattering  conditions  prevail, and  that  the  quadratic 
approximation  to  the  phase  structure-function  is  valid. 

This  approximation  was  used  by  Sreenivasiah  et  al . ,  (1976)  and 
by  Sreenivasiah  and  Ishimaru  (1979)  for  the  cases  of  plane  wave  and  beam 
wave  propagation  in  homogeneous  turbulence.  More  recently  the  two-posi¬ 
tion,  two-frequency  mutual  coherence  function  was  obtained  for  spherical 
wave  propagation  using  the  extended  Huygens -Fresnel  principle  (Fante, 
1981).  Although  the  quadratic  structure-function  approximation  can  some¬ 
times  lead  to  difficulties  (Wandzura,  1980)  it  is  appropriate  for  the 
two-frequency  mutual  coherence  function  but  not  for  calculation  of  higher 
moments  of  the  field  (Fante,  1980).  Fante  (1981)  discusses  the  accuracy 
of  the  quadratic  structure-function  approximation  for  the  case  of  atmo¬ 
spheric  turbulence  with  a  Kolmogordv  power  spectrum.  He  has  found  that 
the  accuracy  is  a  function  of  the  irregularity  power  spectrum  and  of  the 
strength  of  the  turbulence  (Private  Communication,  1982),  with  accuracy 
increasing  for  stronger  scattering; 

This  aspect  of  the  work  here  Is  a  generalization  of  an  earlier 
calculation  (Knepp,  1983(b))  valid  only  foh  isotropic  turbulence.  As  in 
the  former  study,  the  results  here  are  special i zed,  to  the  case  of  a  thin 
phase-screen  approximation  to  the  thick  scattering  medium.  With  this 
simplification,  analytic  results  are  obtained  for  the  received  impulse 
response  function  to  a  transmitted  power  delta  function.  Then  results  may 


easily  be  determined  for  the  mean  time  delay  and  time  delay  jitter  for 
strong,  anisotropic  turbulence  in  the  thin  phase-screen  approximation. 
Results  for  these  quantities  in  this  approximation  have  previously  been 
shown  (Knepp,  1983(b))  to  closely  approximate  the  results  for  a  thick 
scattering  layer  for  the  case  of  isotropic  turbulence. 

The  effects  of  Gaussian  antennas  oh  measurements  of  received 
power,  decorrelation  time  (or  distance),  mean  time  delay,  time  delay 
jitter  and  coherence  bandwidth  are  determined.  It  is  shown  that  aperture 
averaging  can  reduce  observed  signal  power,  increase  observed  decorrela¬ 
tion  time  and  be  a  significant  factor  in  reduced  observed  time-of-arri val 
jitter  at  the  antenna  output. 


SECTIv‘  2 
FORMULATION 


In  this  section  the  solution  for  the  two-position,  two-frequency 
mutual  coherence  functioh  is  presented  for  transmitter  and  receiver 
located  on  opposite  sides  of  a  thick  layer  of  anisotropic  electron  density 
irregularities.  In  the  next  section  this  result  is  utilized  to  determine 
the  effect  of  transmitting  and  receiving  antennas  on  measurements  of  the 
received  signal. 

Consider  a  monochromatic  spherical  wave  ECp,z,u>,t)  which  origi¬ 
nates  from  a  transmitter  located  at  (0,0, -zt)  and  propagates  in  free 
space  in  the  positive  z  direction  where  it  is  incident,  on  an  ionization 
irregularity  layer  which  extends  from  0  _<  z  <_  L.  After  emerging  from  the 
layer  at  z  *  L,  the  wave  then  propagates  in  free  space  to  a  receiver 
located  at  (0,0,zr).  This  thick  layer  geometry  is. shown  in  Figure  1. 

As  the  wave  propagates,  its  phase  substantially  behaves  as  (-i  <k>z+1u>t)  so 
write 

E(?,z,u),t)  «  U{p,z,w)  exp {i (wt-/<k(z ')>dz ’ ) }  (1) 

where  <k(z)>  is  the  mean  v^ave  number  given  by 


M  1/2  „  „  1/2 
<k(z)>  «  -  (1  -  Nfi/nc)  •  k  (1  -  k2/k2) 
c 


where  c  is  the  speed  of  light  in  a  vacuum,  Ne  is  the  mean  ionization 
density,  nc  is  the  critical  electron  density  and  is  related'  to  the  clas¬ 
sical  electron  radius  r  by  n  »  */(X2r  ),  { r  »  2.82*10- 15  m). 

e  c  e  e 


transmitter 
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Figure  1.  Propagation  of  signals  through  a  disturbed 
transionospherle  channel. 


In  the  case  that  the  transmitted  waveform  is  no  longer  a  mono¬ 
chromatic  wave,  hut  can  be  expressed  as  a  waveform  modulated  on-  a  carrier, 
the  two-position,  two-frequency  MCF  r  is  of  interest,  r  is  important 
for  the  calculation  of  pulse  propagation  in  a  random  medium  and  it  serves 
as  a  basis  from  which  to  calculate  the  important  power  impulse  response 
function  and  its  moments.  Under  the  Harkov  approximation,-  f  is  found  to 
satisfy  (Yeh  and  Liu,  1977;  Knepp,  1983(b)). 
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where  ?s  and  Vj  are  the  gradient  operators  in  the  sum  and  difference 
coordinate  system.  Here  r  =  <U(xi,yi,z,u)i)U*(x2,y2»z»t02)s  which  is 
written  as  rfc.n.z.wd)  after  the  sum  and  difference  transformation.  The 
standard  transformations 

x  =  (x1+x2)/2  c  =  xL-x2 

Y  =  (y,+y2)/2  n  =  y  1~y2 

k  =  (k  +k  )/2  k  «  k  -k 

S  12  d  1  2 

are  used  to  obtain  Equation  3.  In  addition,  it  has  been  assumed  that  the 
frequencies  of  interest  are  much  greater  than  the  plasma  frequency.  The 
quantity  k$  is  the  wavenumber  at  the  carrier  frequency,  k$  *  wq/c. 

The  function  A("p)  is  the  integral  of  the  autocorrelation  func¬ 
tion  of  electron  density  fluctuations,  B^,  in  the  direction  of  propa¬ 
gation 

.  oo 

'  A(p1-p2)  *  /  B^-p^z'Jdz'  (4) 

-oo 

where  "Pi-'P2  *  (x  i**2»y  i_y 2)  and  5  *  so  that 

m  _ 

A(pi-p2)  ■  /  /  exp[iKj_*(prp2)]  ❖?{K|_,K2*0)d2kj_  (5) 

■fll 

where  is  the  power  spectrum  of  electron  density  fluctuations.  Equa¬ 
tions  4  and  5  depend  upon  the  validity  of  the  Markov  approximation  where- 
it  is  assumed  that  the  electron  density  fluctuations  are  delta-correlated 
in  the  direction  of  propagation  (Fante,  1975).  That  is  B^(pi-P2,Zi*Z2) 

*  ^i-02)6(Zi*Zz)  . 


The  unknown  two-frequency  MCF  may  be  written  as  r  =  F^g  (Knepp, 
1983(b))  where  r0  is  the  exact  free  space  solution  in  the  parabolic 
approximation.  Substitution  of  r  =  FiT0  into  Equation  3  and  neglect  of 
near-zone  terms  yields  a  differential  equation  with  terms  containing  the 
factors  3 rx / 3X  and  3T j / 3Y.  Since  the  boundary  condition  at  z  =  0, 
r1(c,n,z=0,a)^)  =  1,  is  independent  of  X  and  Y,  the  derivatives  with 
respect  to  these  quantities  may  be  neglected.  The  change  of  variables  z* 

=  z  +  zt,  9  =  S/z‘,  $  a  n/z',  and  the  additional  substitution  ri  =  r2r3 
where 

r3  -  exp  -  *  A(0)(z‘--z  )k4  (L  -  L)2  (6) 

8  t  pyki  k2y 

and  the  replacements 


enable  one  to  write  the  equation  for  r2  as 


ar2 

3v 


L  ^  r2  -  (u2+A2e2)v2r2  *  o 

v2  \3ti2  3 e2/ 


(7a) 

(7b) 

(8a) 

(8b) 

(8c) 


(9) 


where  A(c,n)  has  been  expanded  in  a  truncated  Tayior  series 

A(c,n)  =  A0  +  A2c2  +  A2A2n2  (10) 

(see  appendix).  This  is  the  quadratic  phase  structure  function  approxima 
tion.  The  difference  wave  number  has  also  been  neglected  with 
respect  to  ks.  The  effect  of  this  assumption  is  to  restrict  the  valid¬ 
ity  of  the  solution  to  a  small  range  of  wavelengths  centered  about  ks. 

Thick  Layer  Solution 

An  analytic  solution  of  the  form 

r2  -  f(v)  exp{-g(v)u2-h(y)e2}  (11) 

may  be  substituted  into  Equation  9  to  obtain  three  equations,  the  first 
consisting  of  terms  that  are  factors  of  u4,  the  second  consisting  of 
factors  of  e2,  and  the  third  independent  of  u  and  e.  These  equations  may 
be  solved  exactly  for  r2,  and  the  results  substituted  into  Equation  3  to 
yield 

,  r1(e,4>,z**L+zt ,<*>d)  * 

F  exp {-(A-Bt) 92  -  (A- * )  (12) 


where 


A 


(14) 


-  i(L+zt)  3132 
4 


/  ■  ,  >2  2  2 

(L+z.  )  aia2 

B  = _ l _ 

6 


1  +  3a  xz  tanh  3a  iL 

t  = - - - 

&aiz  +  tanh  3a 


(15) 


(16) 


f 


,1  +  ASa^  tanh  A3a 
A3axzt  +  tanh  Apa^T- 


(17) 


(18) 


In  transcribing  Equation  6  for  r3,  the  relationship  between  the 
coefficient  Ag  and  the  phase  varianpe  given  by  Equation  A-6a  is 
utilized.  Equation  A-5  gives  the  expression  for  the  quantity  A  in  terms 
of  the  geometry  of  the  line-of -sight  relative  to  the  magnetic  field 
direction. 


To  complete  the  solution,  it  is  necessary  to  solve  for  rx  in  the 
region  L  <  z  <  zr  (soe  Figure  1).  Equation  12  serves  as  the  boundary 
condition  at  z1  *  L  +  z^.  Since  the  region  z'  ^  L  *  zt  corresponds  to 
free  space  with  no  Ionization,  Equation  3  is  appropriate  after  the  r0 
substitution  Is  included,  insertion  of  the  z',  9,  and  ♦  substitutions,  and 
deletion  of  the  last  term  involving  the  function  A  that  is  zero  in  free 
space. 
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The  Fourier  transform 


iYfe^.z'.w.)  -■/  /  e1  9 


i  (Kfl9+i< ,  (|>) 


•  X  r1(K9,K^z1,a)£j)  dKQdK# 

may  be  substituted  into  the  suitably  modified  Equation  3  to  obtain  the 
algebraic  equation 


iii  ♦  -  4  ,K9  *  K*)?1  ’  ° 

3z'  2  k  z'2  V  • 

s 


Equation  20  may  be  solved,  and  the  boundary  condition  Equation 
12  applied  at  z'  =  L  +  zt  with  the  result  • 

A  *  ' 

ri(K0,K^,7t+zr,u)d)  =  rl(K9,K*.Uzt.u»d) 


where 


x  exp[-i  Y(K|  ♦  K|)] 


T.l^d  (vl) 

2  k|  (L+zt)(zt+zr) 


The  final  result  may  then  be  obtained  by  taking  the  Fourier  transform  of 
Equation  12  to  obtain  H  OVV  VZr’V  and  then  takin9  the  Averse 
Fourier  transform  of  Equation  21.  The  required  integrals  ane  easily 
found,  and  the  result  may  be  written  as 


I’1(?,n,zt+zr,u)  )  = - ; - 

[  ( 1+i 4 y(A-Bt ) )  (1+i 4y(A-ABt ' ) )  ] 1//2 

C2(A-Bt)/(zt+zr)2  n2(A-ABt')/(zt+zr)2 

l+i4Y(A-Bt)  l+i4Y(A-ABt * ) 

Equation  23  is  the  result  for  n  after  propagation  through  a  thick  layer 
characterized  by  anisotropic  electron  density  irregularities.  The  full 
solution  for  the  two-position,  two-frequency  MCE  is  obtained  by  multi¬ 
plication  by  r0,  the  free  space  MCF.  Since  r0  is  not  affected  by  the 
random  layer,  it  may  be  ignored  here. 

Thin  Phase-Screen  Approximation 

Much  simplification  is  possible  if  the  thick  scattering  layer  is 
replaced  by  an  equivalent  thin  phase-screen  of  infinitesimal  thickness  and 
the  same  overall  phase  variance.  To  the. first  order,  in  the  thin  phase- 
screen  approximation 

i  cosh  AaaiL  3  1 

sinh  AaaiL  *  Ao^iL 

tanh  AaaiL  *  Aaa  tL  (24c) 

Utilising  these  approximations  as  well  as  Equations  8a-c  and  Equation  18 
for  a,,  a2,  and  3,  the  two-position,  two-frequency  MCE  is  found  in  the 
thin  phase-screen  limit  as 


(24a) 

(24b) 
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where 


rl(S,i,zt+zr,»d)  - 


exp 


.2  ,  .2 

<J>  d 


2w 


exp 


l 


z2  = 


2tv*A2 


1 

1/2 

■i  aiUw 

"d  A 

(01  M 

<*>'/&/  _ 

&ll 

n2/(*2/A2) 

la- 

1+i  “d 

a)'/ A2 

’W2r)Ao 


Xz_z^a,A2 


(2b) 


(26) 


(27) 


Equations  25-27  are  valid  for  any  anisotropic  power  spectrum  for  which,  in 

the  strong  scattering  limit,  the  phase  structure-function  can  be  expanded 

in  the  form  of  Equation  10.  For  the  anisotropic  Gaussian  power  spectrum 

of  interest  in  this  work  A  /A,  3  -  r2  where  r  is  the  axial  scale 

o  c  o  o 

size  of  the  elongated  irregularities  as  discussed  in  the  appendix.  The 
parameter  to  is  a  measure  of  the  decorrelation  distance  of  the  complex 
electric  field  as  measured  in  the  plane  of  the  receiver;  u»‘  is  the 
coherence  bandwidth  and  will  be  shown  to  be  inversely  proportional  to  the 
time  delay  jitter;  0$  is  given  in  the  appendix  as  the  phase  standard 
deviation  imposed  on  the  wave  by  the  phase-screen. 


SECTION  3 

APERTURE  ANTENNA  EFFECTS 


In  this  Section,  appropriate  expressions  are  derived  to  obtain 
the  effect  of  the  antenna  aperture  on  measurement*  of  the  properties  of 
the  received  signal  during  strong  scintillation  conditions.  The  first 
geometry  considered  is  that  of  a  spherical  wave  that  propagates  through 
the  disturbed  layer  to  a  receiving  aperture  antenna.  Results  for  this 
case  are  then  easily  extended  to  the  case  of  a  transmitting  and  receiving 
antenna  and.  then  to  the  case  of  a  monostatic  radar  geometry. 

Consider  the  geometry  shown  in  Figure  2  where  the  aperture  is 
located  in  the  receiver  plane,  z  *  zr.  If  the  incident  field  in  the 
plane  of  the  antenna  is  UCp,zr>u),  then  the  complex  voltage  envelope  at 
the  antenna  output  may  be  expressed  as  (Price,  Chesnut  and  Burns,  1972) 

v(?  ,z  ,w)  *  /  'J(p',z  ,w)  A*(p’-p  )  dV  (28) 

or  r  r  o 

where  U  is  the  solution  to  the  parabolic  wave  equation  at  location  H' 
in  the  receiver  plane  for  a  monochromatic  signal  of  frequency  «*,  "p  is 

*  r 

the  location  of  the  center  of  the  aperture  antenna  and  Ap  is  the  complex 
antenna  weighting  function. 

Equation  28  gives  the  output  voltage  from  an  aperture  antenna 
pointed  in  the  z-directlon  for  a  transmitted  monochromatic  waveform  of 
radian  frequency  w.  If  the  transmitted  waveform  is  a  pulse  waveform 
modulated  on  a  carrier,  the  transmitted  signal  can  be  expressed  as 

s(t)  ■  Re  {m( t )  exp(1u  i)}  (29) 


\ 
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Figure  2.  Receiving  antenna  aperture  centred  at  origin  of 
coord i nate  system. 


where  <*^  is  the  carrier  radian  frequency  and  m(r)  is  the  modulation 
waveform.  After  passage  through  a  layer  of  irregularities,  the  received 
signal  may  be  written  (Knepp,  1983(a))  in  the  absence  of  an  antenna,  as 


where 


r(p_.zp,T)  =  Re|e(p,zr,T)exp(itt>oT)| 


e(P.2«.T)  *  —  /  M(ic)U(p,2i>,u+uj  )exp(1ut)dw  (31) 

r  2x.«  r  0 


The  quantity  e  is  called  the  complex  envelope  of  the  received  waveform. 
M(u>)  is  the  Fourier  transform  of  the  transmitted  modulation  waveform  and 
is  given  by 


M( m)  ».  /  m(  t)exp(-1  (Dt)dT 


(32) 


-»  V*  -  « .r  w-  .  »«  .  .  . 
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In  order  to  obtain  the  effect  of  the  ant^md  on  the  received 
time-domain  waveform,  the  received  complex  envelope  U  in  Equation  31  is 
replaced  by  the  expression,  for  v  given  by  Equation  28.  Thus  the 
received  time-domain  signal  at  the  antenna  output  may  be  written  as 


e(p  ,2  ,t)  =  —  I  M(oj)U( p1  ,z  ,w+u)  ) 

0  2ir  r  0 

★  _  _  ^ 

x  A^(  p1 -pQ)exp(i  (DT)d  p’doi 


(33) 


In  the  case  that  the  receiving  aperture  is  an  omnidirectional  point 
antenna.  A*  is  a  delta  function  and  the  received  complex  envelope  given 
by  Equation  33  is  identical  to  the  result  given  by  Equation  31. 


V.3. 

Si* 

m. 


.  *< 


r *  a 


For  the  case  that  the  transmitted  signal  power  has  a  delta 
function  behavior  in  delay,  m(T)m*(t)  =  S(t),  the  correlation  function 
of  the  received  power  can  be  expressed  as 


<efp  ,zr,T)e*(p  ,zr,T)>  =  G_(p  -p  *T) 


•  (2*)4  J  S(K, x) jAr (K) j  exp[iK‘(Pi-P2)]d2K 


(34) 


where  S(K,t)  is  the  generalized  power  spectrum  (Knepp  and  Wittwer,  1984) 
defined  as  the  Fourier  transform  of  the  MCF, 


•V 

M-C 

‘.*A 

L 


y-v, 

f.v.  ■ 


S/K.t)  *  ~  /  r(pd,u»d)  exp  (-iK*pd+i«dT)dwQd2Pd  (35) 

8  IT 


and  A  (K)  is  the  Fourier  transform  of  the  weighting  function  of  the 
r 

receiving  aperture 


A  (K)  *  -4-  /  a„(p)  exp(*iK*p)  d2p 
r  4tt^  r 


(36) 
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The  vector  K  is  given  as  k  sin9(cos<t>  i  +  sin<j>  j)  for  the  geometry 
shown  in  Figure  2. 

To  obtain  Equation  34,  use  is  made  of  the  fact  that  the  MCF 
depends  on  only  the  differences  in  position  and  frequency.  A  detailed 
derivation  similar  to  the  above  appears  in  Ishimaru  (1978).  Note  that  the 
dependence  of  the  power  impulse  response  function,  G,  on  the  receiver 
distance  zr  has  been  omitted  for  notational  convenience. 

Equation  34  can  be  written  in  another  very  useful  form  with  the 
substitution 

f(K,o>  )  =  JL  /  r(7d,<-d)expHK-Pjd2P.  (37) 

4tt 

where  r(K,u)d)  is  the  Fourier  transform  of  the  MCF.  If  r  is  used  in 
Equation  34  one  obtains 

G(pi-p2.t)  *  (2»)3  J  r(K,uJd)|Ar(K)|  2exp[iK*(Pl-p2)+ia)dT]dudd2K  (38) 

It  will  be  seen  that  Equation  38  is  preferable  for  the  evaluation  of  the 
effects  of  apertures,  on  mean  time  delay  and' time  delay  jitter. 

For  later  usage,  it  is  noted  that  . the  function  r(K,u»d)  may  be 
obtained  directly  from  Equation  25  as 


n 


\ 


Aperture  Weighting  Function 

Let  the  antenna  possess  a  Gaussian  beam  pattern  with  gain 

function 

G( 0)  »  Gq  exp  (- 92/ 0^)  (40) 

where  the  receiver  beanwidth  9r  is  related  to  the  effective  aperture 
diameter  D  by  0p  =  (X/D)2/(4£n2)  and  X  is  the  wavelength  at  the 
carrier  frequency.  With  this  choice  of  beamwidth  G(X/2D)  *  Gq/2  so  that 
this  Gaussian  aperture  has  approximately  the  same  3  dB  beamwidth  as  a 
uniformly  illuminated  circular  antenna.  The  antenna  gain  pattern  is  the 
square  of  the  transform  of  the  aperture  response  function  defined  earlier. 
Thus 

Ap(K)  *  GV2  exp  (-02/292)  (41) 

The  Fourier  transform  relationship  of  Equation  36  may  be  used  and  the 
small  angle  approximation  invoked  to  obtain  the  aperture  weighting 
function  and  its  Fourier  transform  as 

Ap(?J  *  2wkW/*  exp  (-kzp202/2) 

and  . 

A.r(K)  *  Qi/2  exp  (-Kz/2k292) 

Before  proceeding,  it  is  useful  to  note  that  the  effect  of  an 
aperture  antenna  is  immediately  apparent  from  the  form  of  Equation  38. 

With  a  point  or  omnidirectional  antenna,  the  beamwidth  is  large  In 
Equation  43  and  Ap(1<)  Is  constant.  Therefore,  with  a  point  antenna,  the 
aperture  angular  response  function  is  not  important  in  Equation  38.  It  is 
apparent  that  the  receiving  aperture  acts  to  modify  the  value  of  T  such 
that 

VVd>  * 


(42) 

(43) 


(44) 


where  the  subscript  a  on  r  denotes  that  this  expression  contains  the 
effect  of  the  recei ving  aperture.  The  subscript  r  on  r  on  the  right 
hand  side  of  Equation  44  refers  to  the  case  here  that  the  receiver  is 
located  at  z^.  This  equation  is  an  explicit  expression  of  the  well  known 
result  that  an  aperture  modifies  the  angular  spectrum  of  the  received 
signal. 

Case  of  Transmitting  and  Receiving  Antennas 

Equations  34  and  38  give  the  received  power  impulse  response 
function  for  the  case  of  an  omni di recti  or  1  spherical  wave  transmitter  and 
a  receiving  aperture  antenna  with  aperture  angular  response  function 
Ar(K).  A  transmitting  aperture  antenna  is  easily  included  in  Equations 
34  and  38  since  its  only  effect  is  to  reduce  the  region  of  the 
phase-screen  that  is  illuminated.  This  is  accomplished  by  including 
an  additional  antenna  angular  response  function  A^Kz^/z^)  where  At(K) 
is  the  response  function  of  the  transmitter  with  beamwidth  0  and  where 
the  factor  zr/zt  transforms  the  antenna  angle  utilized  by  the  trans¬ 
mitter  to  the  angle  observed  by  the  receiver.  This  transformation  is 
equivalent  to  the  observation  that  a  ray  transmitted  from  location  -zt 
at  the  angle  0  and  that  scatters  from  the  thin  layer  is  observed  by  the 
receiver  at  angle  9'  *  z^fir.  For  this  case  one  may  write 

<45> 

where  the  subscript  tr  refers  to  the  fact  that  this  expression  Includes 
the  effects  of  both  transmitting  and  receiving  apertures.  Since  the 
one-way  propagation  path  here  again  utilizes  a  receiver  at  zf,  f  on  the 
right  hand  sine  again  appears  with,  the  subscript  r. 


For  the  later  calculations  of  scattering  loss,  mean  time  delay 
and  time  delay  jitter.  Equation  38  is  evaluated  at  P1-P2  =  0  where 
Equation  45  is  substituted  for  the  factors  r(K,<i>d  )|Ar (K)  | 2.  That  is, 
for  the  one-way  propagation  path  with  both  transmit  and  receive  apertures 

Gtr(0,T)  =  (2ir)3  /  rtr(K,a>d)exp(iuidT)  d2Kd<»)d  (46) 

For  evaluation  of  the  effect  of  an  aperture  antenna  on  the  measured  signal 
decorrelation  distance,  the  necessary  function  is  the  two-position,  single 
frequency  MCF  given  by  the  Fourier  transform  of  Equation  37  as 


r(pd,o)  =  /  rtr(K,0)exp(iK-Pd)d2K 


The  use  of  rtr  includes  the  effect  of  both  transmitting  and  receiving 
apertures. 

Case  Of  Monost.atic  Radar 

Now  consider  the  case  that  the  two-antenna  geometry  above  is 
replaced  by  a  monostatic  radar  geometry.  In  the  radar  case,  the  solution 
is  facilitated  by  considering  two  one-way  propagation  paths.  Let  the 
radar  antenna  be  located  at  (0,0, -?t)  and  also  let  the  antenna  angular 
response  function  be  given  by  A^  with  beamwidth  on  both  trans¬ 
mission  and  reception*  First  consider  the  one-way  propagation  path  from 
the  transmitter  antenna  to  the  target  located  at  (0,0, zr).  For  this 
path,  the  effect  of  the  aperture  antenna  Is  obtained  from  Equation  45  with 
the  omission  of  Ar,  since  the  point  target  accepts  energy  incident  from 
all  directions.  That  is. 
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where  the  subscript  ml  refers  to  the  first  one-way  propagation  path  for 
the  monostatic  radar  geometry. 

Now  consider  the  second  one-way  propagation  path  from  the  target 
back  to  the  receiver.  This  is  the  propagation  geometry  initially 
considered  in  this  report  except  for  the  interchange  of  receiver  and 
transmitter  location.  For  this  one-way  propagation  path,  f  is  replaced 
by  r  appropriate  to  the  one-way  propagation' path  from  zp  to  -z^.. 

Thus 

■  V^-VIV^I2  <49> 

Note  that  the  occurrence  of  here  as  well  as  in  Equation  48 

signifies  the  use  of  the  same  aperture  antenna  (with  beamwidth  0|-)  for 

reception  as  well  as  transmission.  The  only  effect  of  the  reversal  of  the 

direction  of  propagation  is  an  interchange  of  zr  and  zt  in  the 

calculation  of  the  MCF.  This  interchange  has  the  effect  of  replacing  JtQ 

in  Equation  39  with  the  term  l'  =  z.  I  /z  . 

o  t  o  r 

In  the  case  of  the  monostatic  radar  geometry  under  considera¬ 
tion,  the  principle  of  reciprocity  states, that  the  propagating  signal 
takes  the  same  path  over  each  of  the  one-way  propagation  paths.  This  is 
true  provided  that  the  random  medium  is  "frozen"  for  the  duration  of  the 
signal  traversal. 

The  power  impulse  response  function  for  the  signal  received  it 
the  target  due  to  a  transmitted  delta  function  of  power  from  the  radar  is 
given  by  Equation  38  as 

G^O.  t),  «  (2*)3  /  rml(K,a>d)exp{iwdT)d2l(du.d 
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(50) 


If  a  delta  function  of  power  is  transmitted  from  the  target,  the 


power  impulse  response  function  at  the  radar  receiver  is  given  by  a 
similar  expression 

G2(0,t)  =  (2ir) 3  /  rm2(K,a>d)exp(iudT)d2Kdtiid  (51) 

The  net  power  impulse  response  function  after  transmission  of  a  delta 
function  from  the  antenna,  scattering  from  the  target,  and  propagation 
over  both  of  the  one-way  paths  is  given  as  the  convolution  (Ishimaru, 


1978) 

Gm(0,-r)  =  /G1(0,t,)G2(0,T-t*)dt'  (52) 

where  the  subscript  m  refers  to  the  monostatic  radar  geometry.  Using 
the  convolution  theorem,  it  is  easy  to  show  that 

Gm(0,r)  -  (2*)7  /  ?ml(K,u.d)d2K  ^(K1  .«d)d2K*  expO^tJd^  (53) 


This  equation  is  the  equivalent  of  Equation  46  for  the  case  of  the  one-way 
propagation  path. 

The  above  considerations  apply  to  the  calculation  of  the  power 
impulse  response  function  0^(0, t)  for  the  da$e  of  monostatic  radar. 

For  this  geometry  the  signal  decorrelation  distance  is  easily  obtained  as 
a  simplification  of  the  previous  development  for  the  one-way  propagation 
path.  For  a  monostatic  radar,  the  decorrelation  distance  is  obtained  by 
transmission  of  a  signal  to  the  target  and  observation  of  the  scattered 
signal  at  two  spatial  locations  in  the  plane  of  the  receiver.  The  point 
target  acts  as  a  transmitter  of  a  spherical  wave  which  is  then  observed  by 
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the  radar  receiver.  For  a  small  target  and  frozen  irregularities,  there 
is  no  dependence  of  signal  decorrelation  distance  on  the  one-way  propaga¬ 
tion  path  from  the  transmitter  to  the  tarr«t.  The  necessary  two-position, 
single  frequency  MCF  is  given  in  a  form  similar  to  Equation  47  as 

r(pd,0)  =  /  rm2(K,0)exp(iK*Pd)d2K  (54) 

where  rm2  is  given  by  Equation  49  and  includes  the  effect  of  the 
receiving  radar  antenna,.  It  will  be  shown  that  evaluation  of  the  analytic 
form  of  the  expression  rCp^.O)  immediately  yields  the  signal  decorrela¬ 
tion  distance  in  this  strong  scatter  calculation. 
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SECTION  4 
RESULTS 


In  this  section  the  derivations  presented  in  Section  3  are 
utilized  to  obtain  analytic  expressions  for  several  properties  of  the 
received  signal  as  measured  at  the  output  of  the  appropriate  aperture 
antenna. 

Angular  Scattering  Loss 

An  antenna  aperture  acts  to  coherently  collect  the  energy  inci¬ 
dent  upon  the  antenna  and  to  deliver  it  to  the  receiver.  In  the  transmit¬ 
ting  mode  a  directive  antenna  is  designed  to  transmit  energy  only  over  a 
selected  angular  region.  In  the  receiving  mode  this  same  directive 
antenna  accepts  energy  only  from  a  narrow  range  of  angles.  Thus,  relative 
to  an  omnidirectional  point  antenna,  a  directive  antenna  will  experience 
what  is  referred  to  as  angular  scattering  loss  when  the  signal  at  the 
aperture  exhibits  scintillation.  This  angular  scattering  Toss  arises  from 
angle-of -arrival  jitter  present  in  the  incident  wavefront  that  may  cause 
energy  to  arrive  at  the  antenna  propagating  at  angles  greater  than  those 
accepted, by  the  receiving  aperture. 

A  different  but  equivalent  way  to  view  the  effect  of  a  receiving 
antenna  aperture  is  as  a  coherent  integrator  of  the  signal  arriving  on  the 
aperture  face.  If  the  antenna  is  pointing  towards  the  source  of  an  undis¬ 
turbed  plane  wave,  then  the  antenna  output  is  maximum.  If  there  are  fluc¬ 
tuations  in  the  signal  phase  across  the  aperture  or  the  incident  wavefront 
Is  tilted  relative  to  the  antenna  bores ight  direction,  the  coherent  Inte¬ 
grator  output  is  decreased. 


It  is  convenient  to  compute  the  angular  scattering  loss  as  the 
ratio  of  the  total  power  received  with  an  antenna  (or  antennas)  to  that 
received  with  an  omnidirectional  antenna  that  experiences  no  loss.  In 
both  cases  the  transmitted  signal  power  is  taken  as  a  delta  function  in 
delay.  The  angular  scattering  loss  may  be  written  as 

lim 

9  ,9  f  G(0,t)dT 

loss  =  r -  (55) 

/  G(0,i)dT 

where  the  power  impulse  response  function  G  is  given  by  Equation  46  for 
the  one-way  path  with  both  transmit  and  receive  antennas  and  by  Equation 
52  for  the  monostatic  radar  case. 
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The  integral  of  Equation  46  with  respect  to  t  is  easily 
evaluated  using  the  identity  /  exp(i =  2u6(u>d).  This  leaves  a  two 
dimensional  integral  of  r(K,0)  over  angle  which  involves  easily  found 
integrals  of  Gaussian  functions.  The  result  may  be  written  as 

♦  i !k* 

•? 

where  oqx  and  <J0y  are  the  standard  deviations  of  the  theoretical 
angle-of-arri val  fluctuations  in  the  x-  and  y-directions  over  the  propaga¬ 
tion  path  from  -zt  to  zr.  The  quantities  <?0x'  and  a9y*  are 
the  same  standard  deviations  calculated  for  the  one-way  propagation  path 
from  zP  to  -zt  (see  appendix).  The  quantities  and  9r 

are  the  3  dB  beamwidths  of  the  transmitting  and  receiving  antennas, 
respectively. 
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The  expression  for  angular  scattering  loss  given  by  Equation  56 
contains  the  effects  of  variations  of  the  properties  of  the  propagation 
medium  and  of  the  geometry  of  the  propagation  path.  To  show  this  result 
pictorial ly,  let  the  antenna  beamwidths  be  identical  and  let  the  equiva¬ 
lent  phase-screen  be  located  midway  between  transmitter  and  receiver.  In 
this  case  To  simplify  matters  further,  write  the  antenna  beam- 

width  in  terms  of  the  diameter  as  discussed  earlier.  Write  the  angle-of- 
arrival  jitter  as  a  function  of  the  signal  decorrelation  distance  using 
the  expressions  given  in  the  appendix.  These  substitutions  give 

2 a2  /02  «  0.28  D2/*2  and  2 o?  /e2  =  0.28  D2/(£2/A2).  In  order  to 
0x  o  o  ay  o  o 

separate  the  geometric  effects  of  inclination  angle  variations  from 
aperture  averaging  effects,  define  fcp  as  the  decorrelation  distance 
for  propagation  parallel  to  the  magnetic  field  line.  The  quantity  Zp 
is  invariant  with  respect  to  changes  in  the  inclination  angle  and  is  given 
from  Equation  26  and  A-6b  by  Zp  =  /A*q.  The  resu*ts  reraa<n  a  function  of 
the  axial  ratio  q  which  is  here  taken  as  15  as  suggested  by  Wittwer 
(1979). 

Figures  3(a-b )  show  the  angular  scattering  loss  in  decibels  as  a 
function  of  the  relative  antenna  size,  ttytp  and  of  the  inclination 
angle.  In  Figure  3(a)  the  inclination  angle  between  the  magnet i  :  fiel<J 
and  the  direction  of  propagation  is  shown  for  values  of  0°,  15®.  30°,  45® 
and  90°..  In  Figure  3(b)  the  relative  antenna  size  takes. the  values  of  1, 
3,  10,  30  and  100.  It  is  seen  in  both  figures  that  only  when  th»  antenna 
diameter  Is  large  with  respect  to  the  decorrelation  distance  is  the 
angular  scattering  loss  significant.  Small  inclination  angles  ctuse 
increased  values  of  the  observed  phase  standard  deviation  add  thus 
give  increased  angular  scattering  and  increased  scattering  loss  is  shown. 
An  increase  in  the  inclination  angle  causes  a  decrease  in  and  an 
effective  increase  in  the  decorrelation  distance  t©;  therefore,  the 
amount  of  angular  scattering  loss  decreases. 
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INCLINATION  ANGLE  (DECREES) 

(b) 

Angular  scattering  loss  for  one-way  path  with  transmit  and 
receive  aperture  antennas;  a)  versus  relative  antenna  size 
b)  versus  Inclination  angle. 


For  the, case  of  a  monostatic  radar,  the  power  impulse  response 
function  is  given  by  Equation  52  as  the  convolution  of  the  power  impulse 
response  functions  of  each  of  the  one-way  propagation  paths.  A  simple 
change  of  variables  is  utilized  to  reduce  the  integral  /G^O.tJdT  to 
the  product  of  two  integrals  of  Gi  and  G2.  These  integrals  are  easily 
performed  using  the  delta  function  identity  discussed  previously  with 
resulting  angular  scattering  1ors 


where  0t  is  the  radar  antenna  3  dB  beamwidth.  A  comparison  of  this 
result  to  Equation  56  shows  that  the,  addition  of  a  second  one-way  propaga¬ 
tion  path  essentially  squares  the  scattering  loss. 


<  -  S 


Figures  4(a-b)  show  the  angular  scattering  loss  for  the  case  of 
a  monostatic  radar  geometry  as  a  function'  of  the  inclination  angle  and  of 
the  relative  antenna  size  D/tp.  Here  fcp  is  the  decorrelation 
distance  observed  for  propagation  parallel  to  the  magnetic  field  along  the, 
path  from  (0,0,zr)  to  (0,0,-zt).  The  overall  behavior  of  scattering 
loss  is  quite  similar  to  that  exhibited  in  Figure  3.  The  major  difference 
observed  is  the  increased  loss  caused  by  the  additional  one-way  propaga¬ 
tion  path  from  the  target  back  to  the  radar. 

Signal  Oecorrelatlon  Distance 

As  another  aspect  of  the  aperture  aveiaging  effect,  an  aperture 
can  act  to  increase  the  measured  decorrelation  distance  at  the  antenna 
output.  This  observed  increase  is  caused  by  the  averaging  effect  of  the 
coherent  integration  that  smooths  the  more  rapid  fluctuations.  To  obtain 
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the  signal  decorrelation  distance  for  the  one-way  propagation  path  from 
f0,0,-3rt )  to  (G,0,zr),  it  is  necessary  to  perform  the  integration 
specified  in  Equation  47.  Fortunately,  f  given  by  Equation  45  is  the 
product  of  simple  Gaussian  functions.  The  integral  is  available  in 
standard  tables  and  the  resulting  expression  contains  the  factor 

/*ox  "  n  /* gy)  where  the  decorrelation  distances  in  the  x-  and 
y-directions  (see  Figure  A-l)  are 
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If  small  apertures  are  used  so  that  the  beamwidths  are  large 

with  respect  to  the  theoretical  angle-of -arrival  standard  deviations,  the 

parenthetical  expressions  above  are  unity  leaving  l  and  l  in 

3  ox  oy 

agreement  with  the  coefficients  of  c  and  n  given  in  Equation  25  with 
“•d  set  to  zero. 

Figures  5  and  6  show  the  relative  signal  decorrelation  distances 
anrf  *oy^*o^  x~  and  y-directions  for  the  case  of 

identical  transmit  and  receive  antennas  Separated  by  a  centrally  located 
scattering  layer.  The  quantities  t0  and  t0/A  are  the  values  of 
decorrelation  distance  in  the  x-  and  y-directions,  respectively,  that  are 
observed  with  omnidirectional  antennas.  In  the  figures  the  results  are 
again  shown  as  a  function  of  the  ratio  of  antenna  diameter  to  the 
decorrelation  distance  for  propagation  parallel  to  the  magnetic  field 
direction. 
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Figure  6.  Relative  signal  decorrelation  distance  In  the  y-dlrectl 
for  one-way  path  with  transmit  and  receive  aperture  . 
antennas;  a)  versus  relative  antenna  size;  b)  versus 
Inclination  angle. 


Figures  5(a-b)  show  the  relative  decorrelation  distance  in  the 
x-di  recti  on  as  functions  of  D/tp  and  inclination  angle  i|>.  Apertures 
with  diameters  less  than  the  decorrelation  distance  do  not  affect  measure¬ 
ments  of  the  signal  decorrelation  distance.  It  is  evident  that  apertures 
large  with  respect  to  ip  can  greatly  increase  the  measured  decorrela¬ 
tion  distance  (relative  to  an  omnidirectional  antenna)  for  sufficiently 
strong  angular  scattering. 


It  is  immediately  evident  that  the  decorrelation  distance  is  not 

reciprocal.  With  omnidirectional  antennas  the  decorrelation  distances  in 

the  two  orthogonal  directions,  and  A  ,  are  dependent  on  the  path 

geometry  as  given  by  Equation  26.  These  values  are  not  reciprocal; 

i  and  i  are  measures  of  the  average  distance  between  fades  in  the 
ox  oy 

receiver  plane  and  consequently  depend  on  geometry. 


In  the  case  of  the  monostatic  radar,  the  signal  decorrelation 
distances  may  be  obtained  directly  from  Equations  58  and  59  by  using  the 
parameters  appropriate  to  only  the  second  one-way  propagation  path  with 
.transmitter  (target)  at  (0,0, zr).  The  results  are 
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where  i*  and  A'/A  are  the  signal  decorrelation  distances  in  the  x-  and 
oo 

y-  directions  that  correspond  to  the  one-way  propagation  path  from  a 
source  in  the  ir  plane  to  a  receiver  in  the  -*t  plane.  Figures 
7-8  show  Aox  and  t0 y  for  a  monostatic  radar  as  a  function  of 
normalized  antenna  diameter  and  the  angle  between  the  magnetic  field  and 
the  direction  of  propagation.  The  results  are  similar  to  those  shown  in 
Figures  5-6. 
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Mean  Time  Delay  and  Tine  Delay  Jitter 

In  the  same  manner  that  an  aperture  antenna  with  a  small  beam- 
width  neglects  or  averages  out  energy  incident  at  off-boresight  angles  to 
reduce  the  received  power,  it  also  acts  to  reduce  the  measured  time  delay 
and  time  delay  jitter.  This  reduction  occurs  because  the  energy  arriving 
from  directions  away  from  boresight  typically  travels  over  longer  paths 
than  the  more  direct  signal.  These  longer  paths  require  a  longer  propaga¬ 
tion  time  for  that  portion  of  the  received  signal  and  hence  contribute  to 
increased  <T>  and  <*t  values.  If  this  energy  is  neglected  by  an 
aperture  antenna,  then  the  signal  at  the  output  will  be  characterized  by 
smaller  <t>  and  Oj  than  that  measured  by  an  omnidirectional  (point) 
antenna.  •  ■ 

Now  define  the  mean  time  delay,  <t>,  as 
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I  G(0,T)dt 
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and  define  the  time  delay  jitter,  aT,  by  the  second  moment 

„2  .  /  G(0,T)T»dT  _  <t>2  (H) 

!  G(0,t)dT 

where  G(0,t)  is  the  power  impulse  response  measured  at  the  jutput  of  the 
receiving  antenna  due  to  an.  impulse  of  power  originating  from  the  trans¬ 
mitter. 

To  compute  <t>  and  oT  for  the  one-way  propagation  path 
with  transmit  and  receive  aperture  antennas,  it  is  convenient  to  utilize 
Equation  46  for  the  power  impulse  response  function.  The  integrals  of  r 
and  t2  give  first  and  second  derivatives  of  delta  functions  which  yield 
first  and  second  derivatives  of  (*,1*^).  derivatives  are  . 

evaluated  at  zero  ^  as  a  result  of  the  t  integration.  The  remain¬ 
ing  integration  with  respect  to  1C  is  easily  performed.  The  results  may 
be  wri tten  as 


> 
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(64) 


(65) 


Figure  9  shows  the  effect  of  aperture  antennas  on  the  observed 
mean  time  delay  by  showing  the  wean  time  delay  normalized  by  the  value 
observed  with  omnidirectional  antennas,  (1+A2 )/ (2u>* ) .  In  the  figure, 
transmit  and  receive  apertures  are  identical  and  the  equivalent  phase- 
screen  is  assumed  to  be  located  midway  between  the  antennas.  Figure  9(a) 
shows  the  relative  mean  time  delay  as  a  function  of  D/fcp  for  values  of 
the  inclination  angle  of  0°,  15°,  30°,  45°  and  90°.  It  is  seen  that  large 
apertures  can  have  a  significant  effect  in  reducing  the  observed  mean  time 
delay,  particularly  when  the  direction  of  propagation  is  close  to  the  . 
magnetic  field  direction.  In  Figure  9(b)  the  abscissa  and  the  parametric 
quantity  are  interchanged  relative  to  Figure  9(a). 

Figures  10(a-b)  show  the  effect  of  aperture  antennas  on  the 
observed  time  delay  jitter  oT,  Here  is  normalized  to  its  value 
for  omnidirectional  transmit  and  receive  antennas  [  (l+A4)/(2w’2)]  l/,Z. 

It  is  further  assumed  that  the  ratio  o./u  •  is  small,  as  is  always  the 

9  0 

case  for  GHz  and  higher  frequencies.  At  9  *  0®  the  antenna  aperture  has 
identical  effects  on  and  <t>.  At  other  values  of  the  inclination 
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Figure  9; 
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*»)  versus  Inclination  angle.  "M«t1ve  antenna  size; 


RELATIVE  TIME  DELAY  JITTER  RELATIVE  TIME  DELAY  JITTER 


(b) 

Figure  10.  Relative  time  delay  jitter  for  one-nay  path  with  transmit 
and  receive  aperture  antennas;  a)  versus  relative  antenna 
size;  b)  versus  Inclination  angle. 
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Figures  11-12  show  the  values  of  mean  time  delay  and  time  delay 
jitter  for  an  aperture  antenna  normalized  to  their  values  for  an  omni¬ 
directional  antenna.  Again  the  factor  c^/n'g  is  assumed  small  for 
this  monostatic  radar  case.  These  normalized  results  are  very  similar  in 
form  to  those  presented  in  Figures  9-10  and  again  show  the  dramatic  effect 
of  a  large  aperture  antenna  that  operates  in  strong  scattering  conditions. 

Coherence  Bandwidth 

The  coherence  bandwidth  is  appropriately  chosen  as  the  inverse 
of  the  time  delay  jitter  exhibited  by  the  received  waveform  due  to  a 
transmitted  power  that  is  a  delta  function  in  delay  (Knepp,  1983(b)).  In 
the  cases  considered  here  where  o^/cuo  is  small,  the  effects  of 
antenna  apertures  on  the  measured  coherence  bandwidth  are  easily  obtained 
from  Equations  65  and  69.  The  results  for  the  one-way  propagation  path 
with  transmit  and  receive  antennas  and  for  the  monostatic  radar  case  are 
given  as  the  reciprocal  of  the  values  shown  in  Figures  10  and  12, 
respectively. 
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Figure  11. 
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Relative  mean  tine  delay  for  monostatic  radar;  a)  versus 
relative  antenna  size;  b)  versus  Inclination  angle. 
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APPENDIX 


IONIZATION  IRREGULARITY  DESCRIPTION 

The  ionization  irregularities  are  assumed  to  be  elongated  along 
the  direction  of  the  earth's  magnetic  field  as  shown  in. Figure  A-l.  The 
electromagnetic  wave  propagates  in  the  negative  z  direction.  The  magnetic 
field  vector  lies  in  the  y-z  plane  at  an  angle  of  i|»  with  respect  to  the 
z  axis.  The  ionization  irregularities  are  assumed  rotationally  symmetric 
with  autocorrelation  function 


Figure  A-l.  A  single  Irregularity  elongated  along  the  magnetic  field 
line  In  the  y-z  plane. 


where  s  is  measured  along  the  magnetic  field  direction  and  r  is 
measured  perpendicular  to  this  direction.  The  quantity  q  is  known  as 
the  axial  ratio  (Briggs  and  Parkin,  1963)  and  ro  is  the  correlation 
distance  or  irregularity  scale  size. 

The  s,  r  coordinate  system  may  be  related  to  the  x,  y,  z 
system  by  the  equations  s2  *  (y  sin  <J»  +  z  cos  *)2  and  r2  =  x2  + 

(y  cos  -  z  sin  4»)2.  Using  the  above  transformation  the  irregularity 
correlation  function  may  be  expressed  in  the  x,  y,  z.  This  expression  can 
easily  be  integrated  according  to  Equation  6  to  obtain  the  function 
A(c,n)  as 


A(t,n)  = 


a2  A  qrf 
e 


<Ne>2(q2sin2'|>+cos2t|')  ^2 


exp  {  - 


rj  r  2(q2sinz'|*cos  *♦) 


(A-2) 


The  above  equation  is  essentially  Equation  10  in  Briggs  and  Parkin 
(1963).  A(c,n)  is  easily  expanded  in  a  Taylor  series  and,  retaining  only 
terms  up  to  the  quadratic,  one  obtains  Equation  10  in  the  text  where 


A,  . - 1 - ■_  (A— 3) 

<N  >2 
e 

-  q2  / w  qA 

A2  *  - £ -  (A-4) 

<N  >2r 
e  o 
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A  = _ 


q2sin2^  +  cos2'!' 


(A-5) 


This  same  formalism  may  also  be  applied  to  irregularity  power  spectra  that 
are  non-Gaussian.  For  the  case  of  a  power-law  PSD  the  coefficients  A0 
and  A2  are  different  than  above  but  behave  in  essentially  the  same 
manner  as  a  function  of  the  outer  scale,  size  as  long  as  the  three-dimen¬ 
sional  in-situ  electron  density  PSD  falls  off  at  least  as  rapidly  as  K-4. 

Phase  Standard  Deviation 

For  an  ionized  medium,  the  phase  standard  deviation  is  given  as 
the  integral  of  the  irregularity  autocorrelation  function  along  the 
direction  of  propagation  (Salpeter,  1967),  The  result  may  be  written  as  . 
(Knepp,  1983(b)) 


-  k^LA  /k2 

4  p  0  0 


(A-6a ) 


/»  (Xre)2qAroLa2 


(A-6b) 


Theoretical  Angle-Of -Arrival  Fluctuation 


Consider  a  plane  wave  traveling  in  the  negative  z-dlrection  and 
Incident  in  the  x-y  plane  at  an  angle  9  from  the  z-axis.  The  electric 
field  is  given  by  E(x,y)  *  Eq  exp  {ik(sin8  x  ♦  cos9  z ) }  where  the 
exp{i<ut)  time  dependence  is  suppressed.  The  angle-of-arri val  measured 
along  the  x-axis  is  computed  as 
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This  computation  gives  sin  9  which,  in  the  small  angle  limit,  is 
approximately  9. 

In  the  case  of  interest  here  the  incident  field  is  given  as  the 
solution  to  the  parabolic  wave  equation  and  the  angle-of -arrival  may  be 
measured  in  the  x-di recti  on  as 

8  =  _i_  ,  '  (A-8) 

x  ikU0  3x 


with  a  similar  expression  for  9^.  Here  the  measurement  depends  on  the 
frequency.  At  a  single  frequency  the  rms  value  of  8X  is  related  to 
the  second  derivative  of  the  two-position  MCF  (Papoulis,  1965,  p.  317) 


l  a2r(c,n,(i)d=o) 

k2  3t2 


C*  o=0 


(A— 9 ) 


With  a  similar  equation  for  «9y.  For  the  case  of  propagation  from  a 
transmitter  in  the  -zt  plane  to  a  receiver  In  the  zr  plane,  oqx 
and  «8y  may  be  obtained  from  derivatives  of  Equation  25  in  the  text  as 
«9X  *  and  «9y  *  •'?/lk(t0/A) ]• 

For  the  case  of  propagation  from  a  transmitter  in  the  zP 
plane  to  a  receiver  in  the  -zt  plane,  the  interchange  of  transmitter 
and  receiver  may  be  taken  into  account  by  substituting  V  for  t  where 
i*  *  z  t  /z  .  This  substitution  gives  the  angle-of -arrival  fluctuation 

otor  nr  nr-  f  1 

over  the  reversed  path  as  (  *  *2/{ktp  and  cr^  t «  t^/A)  ]• 
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